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Abstract 



> 

O. 

CN ■ We calculate the perihelion precession 5 for nearly circular orbits 

■ in a central potential V(r). Differently from other approaches to this 
CO problem, we do not assume that the potential is close to the Newto- 
OO ! nian one. The main idea in the deduction is to apply the underly- 
ing symmetries of the system to show that 5 must be a function of 

r • V"(r)/V'(r), and to use the transformation behaviour of 5 in a ro- 

■ 

■ tating system of reference. This is equivalent to say that the effective 
potential can be written in a one-parameter set of possibilities as sum 
of centrifugal potential and potential of the central force. We get the 
following universal formula valid for V'(r) > 



5(r) = 2vr • 



1 

1 



y/3 + r-V"{r)/V'{r) 

It has to be read as follows: a circular orbit at this value r exists and 
is stable if and only if this S is a well-defined real; and if this is the 
case, then the angular difference from one perihelion to the next one for 
nearly circular orbits at this r is exactly 27r + 5(r). Then we apply this 
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result to examples of recent interest like modified Newtonian gravity 
and linearized fourth-order gravity. 

In the second part of the paper, we generalize this universal formula 
to static spherically symmetric space-times 

ds 2 = _ e 2A(r )(it 2 + e 2»(r) dr 2 + ^2 

for orbits near r it reads 

5 = 2it- 1 

[y/3-2r-X'(r)+r-X"(r)/X'(r) 

and can be applied to a large class of theories. 

For the Schwarzschild black hole with mass parameter m > it 
leads to 



5 = 2tt- 




a surprisingly unknown formula. It represents a strict result and is 
applicable for all values r > 6m and is in good agreement with the 
fact that stable circular orbits exist for r > 6m only. For r»m, one 
can develop in powers of m and gets the well-known approximation 

67rm 
o ~ . 

r 

Keyword(s): Perihelion advance 

1 Introduction 

Adkins and McDonnell [1] "calculate the precession of Keplerian orbits under 
the influence of arbitrary central- force perturbations." For some examples 
including the Yukawa potential they present the result as hyp ergeo metric 
function. For nearly circular orbits, they arrive at the formula for the peri- 
helion precession A9 P , [1], eq. (11) 



where G is the gravitational constant, M the mass of the central body, m 
the mass of the moving test body, L the radius of the orbit, and u — 1/r 
the inverted radial coordinate. The potential V is the perturbation of the 
Newtonian potential, so the total potential is then given by V(r) — GMm/r. 
They mention that this formula eq. (1.1) is "equivalent to a formula for 
the nearly-circular precession that has been used by Dvali, Gruzinov and 
Zaldarriaga [2]." 

In the fourth section of [1], the Yukawa potential is applied in the form 
[1], eq. (31) 

V(r) = -exp(-r/X) A>0. (1.2) 

r 

Using the parameter k = L/X they arrive at [1], eq. (33) 

1XCX 

A*p(k, 0) = ~~GMm^ eXp( ~ K) • (L3) 

In the fifth section of [1], the authors apply the fact, that within this ap- 
proach, the famous general relativistic perihelion advance can be reproduced 
by using the first post-Newtonian correction 

T/ . G 2 M 2 mL 

v ( r ) = -i-r- ( ) 

where c is the light velocity. They arrive at [1], eq. (42) 

A^ P (GR) = —— (1.5) 

C Li 

and they also present limits for the value of the cosmological constant by 
comparing theoretical and measured values of the Mercury perihelion ad- 
vance. 

The authors of [2] investigate those kinds of theories which possess a 
linearized form of the field equation of the type 



(n + f(D))g ij = T i 



and calculate e.g. the anomalous perihelion precession for this kind of theo- 
ries by perturbations around the Newtonian potential. 
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In [3], the perturbations in the cosmological context are calculated for 
several scalar-tensor theories of gravitation, and for the different conformal 
transformations the distinction between the Einstein and the Jordan frames 
have been made. They applied the results also to calculate an effective grav- 
itational constant for measurements within the Solar system. In [4] and 
[5], the authors carefully calculate the possible measurable effects of tensor- 
multi-scalar theories of gravitation, including the secular rate of perihelion 
advance. 

Davies [6] deduces the perihelion precession due to a perturbing central 
force on an elliptic Keplerian orbit via a perturbation with the Runge-Lenz 
vector. He mentions that one can mimic the influence of the outer planets 
to the perihelion shifts of the inner ones by replacing each outer planet by a 
ring of same total mass, so that the effective potential can remain rotationally 
symmetric. 

2 Perihelion precession 

A test mass shall move along a periodic orbit in a central potential V(r). We 
look for the perihelion precession of this orbit. 

Without loss of generality the test mass has unit mass, and the motion 
takes completely place in the equatorial plane. We parametrize this plane 
by (r, </?), denote the time by t and use the dot to abbreviate for d/dt. Then 
the Lagrangian reads 

L = \ (r 2 + rV) - V(r) . (2.1) 

We assume that V(r) is a C 2 -function at all values r which belong to the 
orbit. For the orbit (r(t), ip(t)) we define perihel and aphel via 

r\ = min r(t) and r 2 = max r(t) (2.2) 

resp., where < r\ < r 2 < oo is assumed. 

Let be the change of ip(t) during the change from r(t) from one perihel 
to the next aphel. Due to time-reversal invariance, the same ip is also the 
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change of (p(t) from this aphel to the next perihel. For 2tp = 2ir, the orbit is 
exactly closed after one revolution. So it is adequate to define the perihelion 
precession 5 by 

5 = 2(v? - tt) • (2.3) 

For purely radial oscillations, our definition implies ipo = 0, i.e. 8 = —2tt. 
If purely radial motion is excluded from the consideration, then all values of 
5 with 5 > —2n may appear as perihelion precession. For 5 > we call it 
perihelion advance. 

What happens with 5 when the orbit is continuously deformed ? Example: 

Let 

r((p) = 4 + £cos(y?) + cos(2y?) (2.4) 

with some parameter e < 1. For all values e > we get 5 = 0, but at e = 
we get 5 = — 7r. This example shows that 5 does not always continuously 
depend on the orbits. 

However, in the typical cases, 5 is a continuous function and for a given 
fixed V(r), we have S(ri, r 2 ), i.e., the prescription of perihel and aphel 
uniquely determine the perihelion precession. We now define 

6(r )= lim 5(n,r 2 ). (2.5) 

ri,7-2— >ro 

The expression 5(ro, tq) is formally the perihelion precession of an exact 
circular orbit which does not make any sense. So, what is the interpretation 
of the limit in eq. (2.5)? It is just the perihelion precession of nearly circular 
orbits which should be well-defined for those cases where the related exact 
circular orbit is a stable one. 

It is the purpose of the present paper to deduce several formulas for the 
calculation of 5 and to apply them to modified Newtonian gravity. 

3 Nearly circular orbits 

How can we calculate the perihelion precession <5(r ) for the nearly circular 
orbit at r = r ? As we have a second order equation of motion, it should 
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be a function of r , V(r ), V'(r ) and V"(r ) only, where the dash denotes 
d/dr. An exact circular orbit at this r-value is possible if and only if the 
repelling centrifugal force is compensated by an attractive central force, i.e., if 
V'(r ) > 0. 1 Now we start to simplify the problem: Adding a constant to the 
potential does not alter the orbits, so no dependence on V(r ) should appear. 
Similarly we argue as follows: if we multiply the function V by a positive 
constant, then we can compensate this by multiplying t also by a suitable 
positive constant without changing the orbits, therefore the dependence on 
the potential can only be in the form of an expression like 

which is invariant with respect to a multiplication of V by a positive constant. 
Finally, we know that 5 is dimensionless, and here we need the last possible 
argument, r , to produce a dimensionless quantity from it: we define 

. n ■ V"(ro) 

q -^ur- (31) 

Example: We assume V(r) = — 1/r, then q = — 2 according to eq. (3.1); this 
potential is the exact Newtonian gravitational field, where we know that all 
the bounded orbits are exact ellipses with the center r = being located at 
one of their focal points, so we get 5 = for this case. This motivates the 
definition q = q + 2, i.e., 

.-=^ + »- M> 

Then it holds: 5 must be a function of q. As no other dependencies exist, it 
must be a universal function 5[q\ being valid for all potentials, and 5[0] = 
because for the Newtonian theory, q — 0. 

The next step is to find out the exact form of this universal function. 
A first idea to assume exact linearity in q is not justified, because then the 
restriction 5 > —2n deduced in the previous section would not be realized. 



1 This sentence is included to fix the sign convention and to make clear, that V'(r ) 
may be written in the denominator. 
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To find the exact form of this universal function it suffices to insert a 
non-trivial one-parameter set of examples for which the solution is known. 

To this end we discuss how the perihelion advance changes if if is replaced 
by p = k-ip with an arbitrary positive parameter k, but r remains unchanged. 
We get ipo = k ■ <p and with eq. (2.3) then 

6 = -2tt + k ■ (5 + 2tt) = k ■ 5 + 2vr • (k - 1) . (3.3) 

The set of possible 5- values is restricted by 5 > —2ir, and by eq. (3.3), also 
5 > —2tt. For k — 1 we get, of course, 8 — 5. To find the related g-values we 
need the equation of motion which is deduced in the next section. 



4 The equation of motion 

The Lagrangian eq. (2.1) reads 

L = \ (r 2 + r V) - V{r) . 

The angular momentum M is a conserved quantity: 

dL o M . . 

M = —— — r (p , hence (p — — . (4.1) 
dip ' r l 

Radial motion is already excluded, so M ^ 0. Without loss of generality we 
assume M > 0, otherwise we would change the orientation of the r— ip— plane. 
The energy E is also a conserved quantity: 

E = ^(r 2 + rV) +V(r). (4.2) 

Inserting eq. (4.1) into eq. (4.2) we get 

f 2 M 2 
E = j + V(r) + — . (4.3) 

We derive eq. (4.3) with respect to t and divide by r afterwards. Then we 
get the Newtonian force equation 

M 2 

= r + V\r) - — , (4.4) 
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the term with M 2 represents the centrifugal force. A circular orbit implies 
f = 0, so by eq. (4.4), this is possible at r = r for V'(r ) > only. 
To evaluate stability, we define the effective potential as usual: 

V eS (r) = V(r) + ^ (4.5) 

leading to 

M 2 

KAr) = V{r) - — (4.6) 

and 

w = no + ^-- (4.7) 

A circular orbit at r = r requires V c g(r ) = E due to eqs. (4.3), (4.5) and 
Kff( r o) = due to eqs. (4.4), eq. (4.6). This implies 



M(r ) = yjr% ■ V'(r ) 

and 

E(r ) = V(r ) + ^ • V'(r ) . 

A simple calculation shows that the following four inequalities are all equiv- 
alent to each other: 

^M >0 , ^>>o, v SM >o 

and 

r • V"(r ) + 3 • V'(r ) > . (4.8) 

A perturbation of the circular orbit can be parametrized by slightly changed 
initial conditions, or equivalently by slightly changed values of M and E. 
In a first step we restrict to perturbations which have the same angular 
momentum M and a slightly changed energy E instead of E. So we have to 
solve 

E = T - + Kfr(r) . (4.9) 

To get solutions one needs E > E. Thus the problem is now equivalent to a 
one-dimensional motion in the potential V e R . From eq. (4.9) we get 



r = ±v / 2- ^E-V eS (r). 
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Together with ip = M/r 2 we find 

dip ip M 

dr f r 2 -V2-^E-V cS (r)' 

If the equation E = V e s(r) has two solutions n, r 2 near r with r\ < r < r 2 
we get 

M dr 

In the limit E — > we have ri,r2 — > ro- We need a positive finite value for 
ipo in this limit, and this is possible for V^g(ro) > only, i.e., if inequality 
(4.8) is valid. If this is fulfilled, then V e e has a regular quadratic minimum at 
r = r and the limit value of ipo depends on V^(r ) only, not on any higher 
derivatives of V(ro). This strictly confirms the assumption made above that 
the universal formula for 5 does not depend on derivatives of V higher than 
the second one. 2 

In a second step we should also look for perturbations where M is slightly 
changed to M. However, such perturbations can be, due to dM(r )/dr > 
0, rearranged to be perturbations at a slightly changed circular orbit with 
adequately chosen f instead of r , so this does not lead to new conditions. 

Now let (r(t), y?(t)) be a periodic solution and k > a parameter. We 
define f(t) = r(t) and (p(t) = k ■ ip(t). We look for a potential V(r) such that 
(f(t), <p(t)) becomes a solution. With eq. (4.1) we get 

M = k ■ M (4.10) 

and with eq. (4.3) 

- f 2 ~ M 2 

E = - + V(r) + — . (4.11) 

2 As it represents a key point in the deduction, we give also the idea for a third inde- 
pendent proof of this statement; it is meant as pedagogic remark: if one considers the 
analogous problem of motion in a 4-dimensional pseudo-Riemannian space-time, then the 
circular orbits are represented by such geodesies, and the nearly circular orbits are repre- 
sented by the geodesic deviation equation, which itself has the components of the curvature 
tensor as coefficients, i.e., no more than second derivatives of the potentials appear; and 
our classical problem of motion can be given as an adequate limit of space-times. 
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An additive constant to the energy can be compensated by adding a constant 
to the potential, so we may assume E = E. A comparison of eq. (4.3) with 
eq. (4.10) and eq. (4.11) leads to 

~ , . k 2 -M 2 \ r / ^ M 2 
V{r) + = V(r) + 2^ 

i.e. to 

V(r) = V(r) + {1 -^ 2 M \ (4.12) 

This means: here we used the transformation behaviour of 5 in a rotating 
system of reference, which is equivalent to say, that the effective potential 
can be written in a one-parameter set of possibilities as sum of centrifugal 
potential and potential of the central force, and so the knowledge about 5 
for one element of this set suffices to calculate it for all other elements of this 
one-parameter set. 3 

Example: Let V(r) = — 1/r, i.e., again the Newtonian potential, we 
consider nearly circular orbits at r = 1. Then V'(r) = 1/r 2 and with eq. 
(4.4) we get using f = just M — 1. Inserting this into eq. (4.12) we get 

Now we apply eq. (3.3). For k = 1 we have, of course, 5 = 0, so we get 

5 = 2tt • (k - 1) . (4.14) 

With the help of eq. (3.2) and eq. (4.13) we calculate at r = 1 

V"(l) 1 

All values q with q > — 1 can appear this way. We invert this equation and 
get with the assumption k > the result 

A; = -=L=. (4.15) 



3 This is the key point of the deduction: we give as input only the knowledge of S for 
the Newtonian potential at r — 1, then by this one-parameter set of transformations we 
produce the knowledge about 5 for a one-parameter class of potentials, and this knowledge 
suffices to identify the universal function 5[q] uniquely. 
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Inserting eq. (4.15) into eq. (4.14) and removing the tilde we get for arbitrary 
values q > — 1 

Indeed, this 5 covers all values S > —2tt and it is defined for all values q > — 1. 
Applying eq. (3.2) and the condition V'(r ) > already mentioned in section 
3 this is equivalent to the conditions 

r ■ V"(r Q ) + 3 • V'(r ) > and V'(r o )>0. (4.17) 

Equation (4.16) together with eq. (3.2) defines the universal function we 
had looked for in section 2. It should be emphasized that we did not solve any 
integrals or differential equations to deduce it, we only applied the obvious 
symmetries of the system and the knowledge about the absence of perihelion 
precession in Newtonian gravity. 

It is still unclear what conditions for the potential V(r) have to be met 
that a periodic orbit with prescribed values of perihel r\ and aphel r 2 exists. 
To this end let us fix a C 2 -function V(r) and values n, r 2 with < r\ < r 2 . 
Both at r 1 and r 2 we have r = 0, so we get from eq. (4.3) 

M 2 M 2 
V (n) + ^ = V(r 2 ) + w (4.18) 

So one needs the finite version of the condition V'(ri) > 0: 

AV = V(r 2 )-V(n)>0. (4.19) 

By the way, the purely radial oscillations which are excluded here, appear as 
M = in eq. (4.18) and require AV = 0. 

Inserting eq. (4.19) into eq. (4.18) and solving for M we get 



M = r 1 r 2] j^ ! . (4.20) 



Inserting eq. (4.20) into eq. (4.3) at r = r\ we get 



E= W^\-W^\ (4.21) 
ri-r{ 
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That the motion between r\ and r 2 is always possible requires f 7^ in this 
whole interval, 4 so we have to fulfil, see eq. (4.3) 

M 2 

E > V(r) + — for n < r < r 2 . (4.22) 
Inserting eqs. (4.20) and (4.21) into this inequality we get 

• (r-2 — r{) 

representing the finite version of the first of the conditions (4.17); one can 
prove this statement by inserting r 2 = r\ + e into eq. (4.23) and applying 
the limit e — > afterwards. 

Finally it should be mentioned that in some limiting cases, also equality 
instead of a <-relation could lead to some solutions; however, in those cases 
either V(r) fails to be a C 2 -function or the test mass would need infinite time 
to reach the limit, however this fails to represent a periodic motion: and both 
is excluded from our considerations. 



5 Nearly circular orbits - second round 

Now we are ready to formulate the result: Let V(r) be a C 2 -function and 
let r > be a fixed value of the radial coordinate. Then an exact circular 
orbit at this r-value is possible if and only if the repelling centrifugal force is 
compensated by an attractive central force, i.e., if V'(r ) > 0, where the dash 
at V denotes d/dr. This orbit represents a stable one in the sense that small 
perturbations of the initial conditions always lead to periodic oscillations 
around r = r , if and only if r • V"(r ) + 3 • V'(r ) > 0. If both inequalities 
are fulfilled, then the perihelion precession 5(r ) of the nearly circular orbits 

4 A priori, an inflexion point might be possible, too, i.e., for instance r > on an 
interval but a single point where r = 0. However, such a behavior is already excluded by 
our assumption that we consider only periodic orbits; a little bit more sophisticated one 
can argue: such a behaviour would a priori allow simultaneously two solutions of the field 
equations, showing that the Cauchy problem would be ill-defined at this point. 
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at r = r is well-defined and can be calculated by use of eqs. (3.2) and (4.16) 



This equation represents an exact result and is not restricted to potentials 
close to the Newtonian one. 

If we rewrite eq. (5.1) in the form 



then imposing the validity of the inequalities (4.17) is equivalent to require 
that eq. (5.2) represents a well-defined real function. This fact can be ex- 
plained as follows: In the deduction of eq. (5.2) we only used symmetry 
arguments and continuous deformations of the orbits, so in the connected 
component of the Newtonian potential V(r) — — 1/r with V'(r) > all regu- 
larity conditions will be met; as V'(r ) = represents a singular point for eq. 
(5.1), it is also clear from that version of the equation why only V'(r) > is 
allowed. 

To ease comparison with the literature it proves useful to work in the 
inverted radial coordinate u = 1/r. We define W(u) = V(r) and a dash at 
W shall denote d/du. Then it holds 



W\u) = -r 2 ■ V\r) , W"(u) = 2r 3 • V'{r) + r 4 • V"(r) . (5.3) 



Because this inversion is a dual transformation we can exchange u with r 
and simultaneously V with W in eq. (5.3) and get 



V'{r) = -u 2 ■ W'{u) , V"{r) = 2u 3 ■ W'{u) + u A ■ W"{u) . (5.4) 



to 




(5.1) 




(5.2) 




(5.5) 




(5.6) 
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In this coordinate the inequalities (4.17) read, see eq. (5.3) 

u ■ W"(uq) - W'(uq) > and W'(uq) < . (5.7) 

Analogously to eq. (5.2) we can now combine eq. (5.6) with inequalities (5.7) 
to get 

, 1 - 1 ] • (5.8) 

^l-u -[\n(-W'(u ))]' J 

If we have the additional condition that \uq ■ W"(uq)\ <C |W'(uo)|, then we 
get from eq. (5.6) the following approximation for 5: 

S(l/u ) — 7r • Uo • W"{u )/\W'{u )\ . (5.9) 

The Newtonian potential for a central mass to > reads V(r) = —mG/r, 
where G is the gravitational constant. This leads to W(u) = —mGu, hence 
W'(u) = —mG and W"(u) = 0. Therefore, eq. (5.9) is especially useful 
if the central potential under consideration is a small perturbation of the 
Newtonian potential. 



5{l/u Q ) = 2tt- 



6 Comparing with Adkins and McDonnell 

In [1], see also the references cited there, the orbital precession due to central- 
force perturbations has been calculated in details, and applications are given; 
especially, their eq. (11) (i.e. our eq. (1.1) above) 

p - GMmL du 2 lu=1/L 

is, after adequate transformation of the notation, almost identical to our eq. 
(5.9). 

Let us check this statement in more details. To this end we now transform 
our formulas to the notation used in [1]. This first means: from now on we 
denote the central mass by M and the mass of the test body by to, and we 
reintroduce gravitational constant G and light velocity c into the formulas, 
even in those cases, where units have been chosen with G = c = 1. Our W(u) 
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in eq. (5.9) has then to be replaced by V(u) — GMmu, where now this V(u) 
is the small perturbation of the Newtonian potential according to [1]. This 
leads to W'{u) = V\u) - GMm and W"{u) = V"{u). So, for the second 
derivative, there is no difference. For the first derivative, we have within the 
used approximation, that V'(u) can be neglected in comparison with GMm. 
So far, the r.h.s. of eq. (5.9) reads 

nu V"(u ) 
GMm 

and evaluating this at uq — 1/L exactly leads to the r.h.s. of eq. (1.1). 
Example: Let V be according to eq. (1.2), i.e. 

Oi 

V(r) = -exp(-r/A) A>0. 

Then we have 

W{u) = -GMmu ■ [1 - /3exp(-l/(Au))] (6.1) 

with a = GMmf3. The perihelion shift according to [1] is with k — L/X, see 
eq. (1.3): 

A6 P = -7r/^ 2 exp(-/t) (6.2) 

an expression which is, as a consequence of the approximation used, com- 
pletely linear in the parameter (5. In fact, the result is valid only in regions, 
where the perturbation is sufficiently small. 

Now we apply our formula eq. (5.6) to the same problem eq. (6.1). The 
factor GMm will cancel out anyhow in eq. (5.6), so we may put GMm = 1, 
i.e. a = (3, already now. We get 

W(u) = -« + /3u-exp(-l/(A«)), W\u) = -1 + (l + -exp(-l/(Au)) . 

(6.3) 

The more conventional form of this potential W appears when one writes it 
in dependence on r = 1/u to get 

W = -- + -■ exp(-r/A) . (6.4) 
r r 
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In the present case, the second inequality (5.7) evaluated at uq — 1/L and 
using k = L/X > reads 

P ■ (1 + k) < exp(/c) (6.5) 

and gives a restriction for /3 > 1 only. In details: let us fix any Ko > 0, then 
we define 

/3 = exp(«o)/(l + kq) 

and then (6.5) is fulfilled for k > Kq only. 

In other words: For every /3 < 1, all positive radius values L appear for 
a circular orbit. For each j3 > 1, there is a positive r such that a circular 
orbit exists for L > r only. 5 

For the second derivative we get from eq. (6.3): 

W"(u) = -^- 3 ex V (-l/(Xu)). (6.6) 
Inserting eqs. (6.3) and (6.6) into eq. (5.6) we get 

S(L) = 2,n ■ ( 1 k2 - l) . (6.7) 

\ V l+«-exp(«)//9 / 

Let us examine eq. (6.7): For very small values \/3\ 1 it is continuous and 
eq. (6.2) is a good approximation to it in correspondence with the fact, that 
here the perturbation to the Newtonian potential is small and therefore, eq. 
(6.2) is applicable. 

Things are quite different for other cases: From eq. (6.4) one can see 
that only for f3 — 1, the potential remains bounded as r — > 0. For this case 
and small values of k < 1 eq. (6.2) gives 5 = —ttk 2 , whereas the exact 
formula eq. (6.7) has in the same case 5 = -2tt(1 -1/V3), a totally different 
behaviour. 

5 Hcre, this example was chosen mainly to present how to apply our formula. However, 
if one wants to interpret the physics behind, one should note that positive values of (3 
would correspond to ghost degrees of freedom if they have an even helicity. But if the 
extra Yukawa force is mediated by a vector field (like the W and Z bosons), then even 
positive values of (3 are allowed. 



16 



In section 3 of [1] it is mentioned that the development of perihelion 
precession in powers of the eccentricity e contains only even powers of e; 
this means that for sufficiently small values of e, where linearization in e is 
justified, our formulas for nearly circular orbits are applicable, too. 

7 Application to fourth-order gravity 

Further examples are as follows: In [7] I wrote without explicit proof, see also 
[8], page 235-236: "Next, let us study the perihelion advance for distorted 
circle-like orbits. Besides the general relativistic perihelion advance, which 
vanishes in the Newtonian limit, we have an additional one of the following 
behaviour: for r — > and r — > oo it vanishes and for r m l/m® and r ps \jmi 
it has local maxima, i.e., resonances." 

This refers to linearized fourth-order gravity, see Stelle [9] for details, 
where the gravitational potential for a point mass m reads 

V(r) = — mr _1 (l + exp(— m r)/3 — 4exp(— m 2 r)/3) . (7.1) 

The perihelion precession of this and similar theories can be calculated by 
inserting this potential V(r) into the equation (5.1). Of course, in the region 
of large values r, the known approximations like (5.9) would serve also, but 
our equation (5.1) will give the correct result also for those r- values, where 
V(r) is far from being close to the Newtonian potential. 

Here in eq. (7.1), m is the mass of the massive spin 0-graviton stemming 
from the i? 2 -term in the Lagrangian, and m 2 is the mass of the massive spin 
2-graviton 6 stemming from the term CijkiC tjkl in the Lagrangian. Both m 
and rri2 are assumed to be positive to exclude the appearance of tachyons, 
but m — > oo and m 2 — > oo represent sensible limits. 

In the case m = m 2 > 0, eq. (7.1) exactly leads to the case (3=1 
discussed in the previous section, the other cases are similar. 

6 This massive spin-2 excitation is a ghost, i.e., carries negative kinetic energy and 
thereby spoils the stability of the model. Therefore, this contribution has not a direct 
phenomenological interpretation. 
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8 Discussion - first part 



At that time paper [7] was first published, in 1986, there this was a purely 
theoretical question. However, recently there is a development to take such 
quadratic gravity theories quite seriously in the sense that their predictions 
can be confronted with observations, see e.g. [10], [11], [12], [13], [14] and the 
references cited there. Also the cosmological solutions of this kind of theories 
have been analyzed in more new details recently, see e.g. [15], [16], [17], [18], 
[19], [20], [21], [22], [23], [24], [25], [26] and [27]. 

Further, it should be mentioned, that for very distorted orbits, [28] and 
[29] give exact results for the perihelion precession for a perturbed Newtonian 
potential. 

Now, let us look for which potentials V(r) the parameter 5(r ) takes 
values according to equation (5.1) which do not depend on r . After some 
calculation we get: up to the inessential transformations mentioned above, 
there is a parameter c > such that 

5 = 2tt(c- 1/2 - 1) (8.1) 

and 

V(r) = lnr for c = 2 , V(r) = — !— ■ r c ~ 2 else. (8.2) 

c — 2 

As expected, one just gets the self-similar functions as solutions to this prob- 
lem. The case c = 1, i.e., 5 = just recovers the Newtonian case. Here the 
bounded orbits are all exact ellipses, the center of symmetry of the potential, 
r = 0, being at one of their focal points. 

For c = 4, eq. (8.2) leads to the harmonic oscillator V = r 2 /2; here 
the bounded orbits are also ellipses, but now, the center of symmetry of 
the potential coincides with the center of the ellipses, therefore, the next 
perihelion is already after one half rotation, i.e. ip = n and 5 = —it in 
accordance with eq. (8.1) for this case. This result once more confirms that 
our result is a strict one also far from the Newtonian potential. 

To prepare for the next part, we now apply units such that light velocity 
c = 1; then it holds: the velocity of the test particle in the exact circular 
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orbit at r = r is less than light velocity if 

r • V'(r ) < 1 . (8.3) 

9 Circular and nearly circular geodesies 

In this second part of the paper we generalize the results of the first part to 
static spherically symmetric space-times, see also [30], [31], [32] and [33] for 
other papers on similar topics. 

Let us now generalize the resulting eq. (5.1) to the analogous situation in 
a 4-dimensional static spherically symmetric space-time. 7 We additionally 
assume that Schwarzschild coordinates are possible, so we consider the metric 

ds 2 = -e 2X dt 2 + e 2 »dr 2 + rW (9.1) 

where dVL 2 is the metric of the standard 2-sphere and A and \i depend on 
r only. We look for time-like geodesies in this space-time (9.1). After suit- 
able rotation of the coordinate system this geodesic remains completely in 
the equatorial plane. Due to the chosen symmetry it holds: Geodesies in 
the equatorial plane of the 4-dimensional space-time (9.1) are exactly the 
geodesies in the 3-dimensional space-time 

ds 2 = -e 2X dt 2 + e 2 ^dr 2 + r 2 d<p 2 . (9.2) 

The coordinates in (9.2) are x l where % — 0, 1, 2, and the geodesic shall be 
parametrized by its natural parameter r: 

x\t) = (t(r), r(r), tp( T )) . (9.3) 

7 A deduction fully analogous to that one from the first part seems not to be easily 
done. The second variant, namely to apply the geodesic deviation equation, also leads to 
unnecessary complicated expressions. As third idea one could try to apply general exact 
solutions of the geodesic equation as found in the text-book literature, e.g. [34], but the 
elliptic integrals appearing there are not easy to handle, therefore we now choose a fourth 
method, namely the direct calculation with nearly circular geodesies. 
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With a dot denoting d/dr we get from (9.2) 

_l = _ e 2^2 + e2/V 2 + r 2^2 (94) 

We may assume t > 0. The components of the geodesic equation read: 

= i'+2A'ir, (9.5) 
= e 2/1 (f + //r 2 ) + \' e 2X i 2 - np 2 (9.6) 

and 

O = ry3 + 2r0. (9.7) 
We define angular momentum M as usual: 

M = r 2 (p. (9.8) 

Purely radial motion shall not be considered, so we have M^O. Without 
loss of generality we may assume M > 0, for otherwise we could reverse the 
orientation of the equatorial plane. 

Due to eq. (9.7), M is a conserved quantity, and we apply this fact to 
simplify eqs. (9.4) and (9.6) to 

_1 = _ e 2A t 2 + e 2fl r 2 + M 2 /r 2 (9.9) 

and 

= e 2 » (f + fi'r 2 ) + \'e 2X i 2 - M 2 /r 3 (9.10) 
resp. Inserting eq. (9.9) into eq. (9.10) we can cancel t to get 

= e *(* + (Y + ^) + Y(l + ^)-£. ("I) 

Next we look for a circular orbit at a fixed value ro of the radial coordinate 
r(r) = r > 0. With eq. (9.11) we get 

z Ad 2 

r • A'(ro) = , where z = —5- > . (9.12) 

2+1 
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This means: a circular time-like geodesic orbit at r = r exists if and only if 
the inequalities 

A'(ro) > and r • A'(r ) < 1 (9.13) 

are fulfilled. They are fully analogous to the second of the inequalities (4.17) 
and inequality (8.3) resp. 
From eq. (9.12) we get 

M 2 = , 4 ' V ^> ■ . (9.14) 
1 - r • A'(r ) 

For M eq. (9.14) the condition dM/dr > is equivalent to 

r • A"(r ) + 3 • A'(r ) > 2r (A'(r )) 2 . (9.15) 

Next, let us define energy E by 

E = e 2A(r) i > . (9.16) 

Due to eq. (9.5), E is a conserved quantity. We can apply this equation to 
remove t from eq. (9.9), leading to 8 

-1 = -E 2 e" 2A + e 2M r 2 + M 2 /r 2 . (9.17) 

For circular orbits at r = r this leads to 

E 2 = e 2X(r )f 1 + Ml) . (9.1 8 ) 



r, 



Inserting eq. (9.14) into eq. (9.18) we get for the energy of the circular orbit 

at r = r 



e A(r ) 



E = (9.19) 
yJl-ro-X'(r ) 

The condition dE/dr > is equivalent to the condition (9.15). 



8 The material in this section is essentially text-book standard, as can be found e.g. in 
[34]; but we presented it here in details to maintain a self-consistent notation. 
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10 Perihelion precession in space-time 



We now prescribe a value r > such that (9.13) and (9.15) are fulfilled. 
The circular orbit at r = tq has angular momentum according to (9.14) and 
energy E according to (9.19). This circular orbit shall now be perturbed, the 
perturbed orbit shall have the same angular momentum M but a different 
energy E ^ E. For the radial coordinate r in dependence on r we make the 
following ansatz 

r = r + e ■ sm(ar) (10.1) 

where a is a positive parameter and e shall be small such that higher powers 
of e may be neglected. We insert this ansatz (10.1) into eq. (9.11) and get 
the following identity 



a = e -,M . ^o-A w (ro)+3-^(ro)-2ro(^(ro)) 2 ^ 

y/l-r - A'(r ) • y/r^ 

It is remarkable that just the inequalities deduced before ensure that a be- 
comes a well-defined positive real. From eqs. (10.1) and (10.2) we get: The 
time from one perihelion to the next is then To defined by 

r = 27r/a. (10.3) 

The perihelion shift 5 is defined as 

<5 = V?(r )-v?(0)-27r, (10.4) 

it measures how much the change in the angular coordinate tp differs from 
27r when the orbit changes from one perihelion to the next one. Clearly, for 
5 = 0, the orbits are exactly closed after one revolution. From eqs. (9.8), 
(10.3) and (10.4) we get 

c- 2f M . 

5 = t -2tt 10.5 

a r 2 



thus leading to the final result 

pM r o) 

5 = 2?r • 



^3-2r -A'(ro)+r -A"(ro)/A'(ro) 



10.6) 
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11 Discussion - second part 



The final formula eq. (10.6) has a structure quite similar to the corresponding 
formula eq. (5.1) from the first part. But a direct change over from one 
of them to the other one is not easily done, so it was really necessary to 
deduce both of them. In eq. (10.6) one can observe, that the spatial metric 
component encoded by the function \i essentially enters the formula but none 
of their derivatives do enter here. This is in contrast to the temporal metric 
component encoded by the function A from which only the first and second 
derivative do enter. 

As a first test, let us insert the Schwarzschild-de Sitter solution into eq. 
(10.6). In units where c = G — 1, we have to insert into eq. (9.1) 



-,2A 



= = 1 _ *1 - ± . 



;n.i) 



r 3 

where m > is the mass of the source and A > the cosmological constant. 
For A = this leads to 



2tt 



1 



- 1 



(11.2) 



— 6m/r 

This is a strict result and is applicable for all values r > 6m. Eq. (11.2) is 
surprisingly unknown up to now. It is in good agreement with the fact that 
stable circular orbits exist for r > 6m only. For r m, one can develop in 
powers of m and gets the well-known approximation 

67rm 



(11.3) 



Example: for ro = 24m, the exact equation (11.2) leads to 5 — 7r-(4/y / 3— 2) = 
55, 7°, whereas the approximation (11.3) leads to 5 = 7r/4 = 45°. 

For A > 0, the condition that a time-like circular orbit exists, is the same 
as for A = 0, namely r > 3m, but the formula for perihelion shift becomes 
a little bit more complicated, 



2tt- 



^Jl- 6m/r Q 



Arf 
6m 



3 + 



9m 
r — 6m 



11.4) 
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but in the case A << 1/m 2 and r not too large one gets the useful approx- 
imation 

67rm 7rArn 
5 w 1 . 

r m 

In a following paper, eq. (10.6) shall be applied also to other spherically 
symmetric metrics. 
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